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Abstract We review recent progress in applying the AdS/CFT correspondence to finite-temperature field theory. 
In particular, we show how the hydrodynamic behavior of field theory is reflected in the low-momentum limit of 
i-C ■ correlation functions computed through a real-time AdS/CFT prescription, which we formulate. We also show how 

the hydrodynamic modes in field theory correspond to the low-lying quasinormal modes of the AdS black p-brane 
metric. We provide a proof of the universality of the viscosity /entropy ratio within a class of theories with gravity 
duals and formulate a viscosity bound conjecture. Possible implications for real systems are mentioned. 
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This review is about the recently emerging connection, through the gauge/gravity correspondence, 
between hydrodynamics and black hole physics. 

The study of quantum field theory at high temperature has a long history. It was first motivated 
by the Big Bang cosmology when it was hoped that early phase transitions might leave some 
imprints on the Universe pQ. One of those phase transitions is the QCD phase transitions (which 
could actually be a crossover) which happened at a temperature around T c ~ 200 MeV, when 
matter turned from a gas of quarks and gluons (the quark-gluon plasma, or QGP) into a gas of 
hadrons. 

An experimental program was designed to create and study the QGP by colliding two heavy 
atomic nuclei. Most recent experiments are conducted at the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory. Although significant circumstantial evidence for the QGP 
was accumulated [2], a theoretical interpretation of most of the experimental data proved difficult, 
because the QGP created at RHIC is far from being a weakly coupled gas of quarks and gluons. 
Indeed, the temperature of the plasma, as inferred from the spectrum of final particles, is only 
approximately 170 MeV, near the confinement scale of QCD. This is deep in the nonperturbative 
regime of QCD, where reliable theoretical tools are lacking. Most notably, the kinetic coefficients 
of the QGP, which enter the hydrodynamic equations (reviewed in Sec. [2]), are not theoretically 
computable at these temperatures. 

The paucity of information about the kinetic coefficients of the QGP in particular and of strongly 
coupled thermal quantum field theories in general is one of the main reasons for our interest 
in their computation in a class of strongly coupled field theories, even though this class does 
not include QCD. The necessary technological tool is the anti-de Sitter-conformal field theory 
(AdS/CFT) correspondence [3jlU[5], discovered in the investigation of D-branes in string theory. 
This correspondence allows one to describe the thermal plasmas in these theories in terms of black 
holes in AdS space. The AdS/CFT correspondence is reviewed in Sec.[3j 

The first calculation of this type, that of the shear viscosity in N = 4 supersymmetric Yang- 
Mills (SYM) theory [6], is followed by the theoretical work to establish the rules of real-time 
finite-temperature AdS/CFT correspondence 0[8]. Applications of these rules to various special 
cases [91 [TU1 [TT1 [12] clearly show that even very exotic field theories, when heated up to finite 
temperature, behave hydrodynamically at large distances and time scales (provided that the number 
of space-time dimensions is 2+1 or higher). This development is reviewed in Sec. [H Moreover, the 
way AdS/CFT works reveals deep connections to properties of black holes in classical gravity. For 
example, the hydrodynamic modes of a thermal medium are mapped, through the correspondence, 
to the low-lying quasi-normal modes of a black-brane metric. It seems that our understanding of the 
connection between hydrodynamics and black hole physics is still incomplete; we may understand 
more about gravity by studying thermal field theories. One idea along this direction is reviewed in 
Sec.0 

From the point of view of heavy-ion (QGP) physics, a particularly interesting finding has been 
the formulation of a conjecture on the lowest possible value of the ratio of viscosity and volume 
density of entropy. This conjecture was motivated by the universality of this ratio in theories with 
gravity duals. This is reviewed in Sec. EJ 
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This review is written primarily for readers with a background in QCD and QGP physics who are 
interested in learning about AdS/CFT correspondence and its applications to finite-temperature 
field theory. Some parts of this review (for example, the section about hydrodynamics) should be 
useful for readers with a string theory or general relativity background who are interested in the 
connection between string theory, gravity, and hydrodynamics. The perspectives here are shaped 
by our personal taste and therefore may appear narrow, but the authors believe that this review 
may serve as the starting point to explore the much richer original literature. 

In this review we use the "mostly plus" metric signature — h ++. 

2 HYDRODYNAMICS 

From the modern perspective, hydrodynamics [13] is best thought of as an effective theory, describ- 
ing the dynamics at large distances and time-scales. Unlike the familiar effective field theories (for 
example, the chiral perturbation theory), it is normally formulated in the language of equations of 
motion instead of an action principle. The reason for this is the presence of dissipation in thermal 
media. 

In the simplest case, the hydrodynamic equations are just the laws of conservation of energy and 
momentum, 

8^ = 0. (1) 

To close the system of equations, we must reduce the number of independent elements of T^ v '. 
This is done through the assumption of local thermal equilibrium: If perturbations have long wave- 
lengths, the state of the system, at a given time, is determined by the temperature as a function 
of coordinates T(x) and the local fluid velocity u^, which is also a function of coordinates u^(x). 
Because u^u^ = — 1, only three components of are independent. The number of hydrodynamic 
variables is four, equal to the number of equations. 

In hydrodynamics we express T^ v through T(x) and u^(x) through the so-called constitutive 
equations. Following the standard procedure of effective field theories, we expand in powers of 
spatial derivatives. To zeroth order, T^ u is given by the familiar formula for ideal fluids, 

= (e + P)u^u v + Pg» v , (2) 

where e is the energy density, and P is the pressure. Normally one would stop at this leading 
order, but qualitatively new effects necessitate going to the next order. Indeed, from Eq. [2] and the 
thermodynamic relations ale = TdS, dP = sdT, and e + P = Ts (s is the entropy per unit volume), 
one finds that entropy is conserved [H] 

d M (sO = 0. (3) 

Thus, to have entropy production, one needs to go to the next order in the derivative expansion. 
At the next order, we write 

= (e + P)u»u u + Pg^ - , (4) 

where a^ u is proportional to derivatives of T{x) and u^{x) and is termed the dissipative part of 
T^ v . To write these terms, let us first fix a point x and go to the local rest frame where u l {x) = 0. 
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In this frame, in principle one can have dissipative corrections to the energy-momentum density 
T ^. However, one recalls that the choice of T and is arbitrary, and thus one can always redefine 
them so that these corrections vanish, cr 00 = a 0t = 0, and so at a point x, 



T UU = £j jUl = q _ ( 5 ) 

The only nonzero elements of the dissipative energy-momentum tensor are <Jij. To the next-to- 
leading order there are extra contributions whose forms are dictated by rotational symmetry: 



r) (diUj + djUi - ^Sijd k u k ^j + (5ijd k u k . (6) 



Going back to the general frame, we can now write the dissipative part of the energy-momentum 
tensor as 



V (d a U(i + d/3U a - ^gapdxu^ + C,g a pd\u x 



(7) 



where P^ u = + u^u u is the projection operator onto the directions perpendicular to u^. 

If the system contains a conserved current, there is an additional hydrodynamic equation related 
to the current conservation, 

d»f = 0. (8) 
The constitutive equation contains two terms: 

f = pu v - DP^d v a , (9) 

where p is the charge density in the fluid rest frame and D is some constant. The first term 
corresponds to convection, the second one to diffusion. In the fluid rest frame, j = —DX7p, which 
is Fick's law of diffusion, with D being the diffusion constant. 



2.1 Kubo's Formula For Viscosity 

As mentioned above, the hydrodynamic equations can be thought of as an effective theory describing 
the dynamics of the system at large lengths and time scales. Therefore one should be able to use 
these equations to extract information about the low-momentum behavior of Green's functions in 
the original theory. 

For example, let us recall how the two-point correlation functions can be extracted. If we couple 
sources J a (x) to a set of (bosonic) operators O a (x), so that the new action is 

S = S + [ J a (x)O a (x), (10) 

J X 

then the source will introduce a perturbation of the system. In particular, the average values 
of O a will differ from the equilibrium values, which we assume to be zero. If J a are small, the 
perturbations are given by the linear response theory as 

(O a (x)) = - [ GK(x-y)J b (y), (11) 
J y 

where G^ b is the retarded Green's function 

iG*(x -y) = 8(x° - y°)([O a (x), O h {y)\) ■ (12) 
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The fact that the linear response is determined by the retarded (and not by any other) Green's 
function is obvious from causality: The source can influence the system only after it has been 
turned on. 

Thus, to determine the correlation functions of T^ v ', we need to couple a weak source to T^ v and 
determine the average value of T^ u after this source is turned on. To find these correlators at low 
momenta, we can use the hydrodynamic theory. So far in our treatment of hydrodynamics we have 
included no source coupled to T^ v . This deficiency can be easily corrected, as the source of the 
energy-momentum tensor is the metric g^ v . One must generalize the hydrodynamic equations to 
curved space-time and from it determine the response of the thermal medium to a weak perturbation 
of the metric. This procedure is rather straightforward and in the interest of space is left as an 
exercise to the reader. 

Here we concentrate on a particular case when the metric perturbation is homogeneous in space 
but time dependent: 

9ij(t,x) = Sij+hij(t), hij<l (13) 
Sdo(t,x) = -1, g 0i (t,x) = 0. (14) 

Moreover, we assume the perturbation to be traceless, hn = 0. Because the perturbation is spatially 
homogeneous, if the fluid moves, it can only move uniformly: u l = u l (t). However, this possibility 
can be ruled out by parity, so the fluid must remain at rest all the time: = (1,0,0,0). We now 
compute the dissipative part of the stress-energy tensor. The generalization of Eq. [7] to curved 
space-time is 



a fJ.v — pfiapv/3 



r](V a up + VpUa) + ( C - -1} ) g aj3 V ■ u 



(15) 
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Substituting u M = (1,0, 0, 0) and g^ v from Eq. [131 we find only contributions to the traceless spatial 
components, and these contributions come entirely from the Christoffel symbols in the covariant 
derivatives. For example, 

a xy = 2nT xy = nd h xy . (16) 

By comparison with the expectation from the linear response theory, this equation means that we 
have found the zero spatial momentum, low-frequency limit of the retarded Green's function of 

rpxy. 

G*, xy (u,0) = Jdtdxe iujt 9(t)([T xy (t,x), T xy (0,0)}} = -ir,u + 0{oj 2 ) (17) 

(modulo contact terms). We have, in essence, derived the Kubo's formula relating the shear viscosity 
and a Green's function: ^ 

77 = - lim -ImGj ' (w,0). (18) 

There is a similar Kubo's relation for the charge diffusion constant D. 



2.2 Hydrodynamic Modes 

If one is interested only in the locations of the poles of the correlators, one can simply look for the 
normal modes of the linearized hydrodynamic equations, that is, solutions that behave as e - iaJ *+ lk - x . 
Owing to dissipation, the frequency o>(k) is complex. For example, the equation of charge diffusion, 

d t p - DV 2 p = 0, (19) 
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corresponds to a pole in the current-current correlator at u = —iDk 2 . 

To find the poles in the correlators between elements of the stress-energy tensor one can, without 
loss of generality, choose the coordinate system so that k is aligned along the x 3 -axis: k = (0, 0, k). 
Then one can distinguish two types of normal modes: 

1. Shear modes correspond to the fluctuations of pairs of components T 0a and T 3a , where a = 
1,2. The constitutive equation is 

T 3a = -nd 3 u a = -^_d 3 T 0a , (20) 

and the equation for T 0a is 

d t T 0a - j^pdiT 0a = . (21) 

That is, it has the form of a diffusion equation for T 0a . Substituting e - tujt + tkxS into the 
equation, one finds the dispersion law 

u = -i-^—k 2 . (22) 

e + P y J 

2. Sound modes are fluctuations of T 00 , T 03 , and T 33 . There are now two conservation equations, 
and by diagonalizing them one finds the dispersion law 

- = ^-Kl" + c )7TF- <23> 



where c s = y/dP/de. This is simply the sound wave, which involves the fluctuation of the en- 
ergy density. It propagates with velocity c s , and its damping is related to a linear combination 
of shear and bulk viscosities. 

In CFTs it is possible to use conformal Ward identities to show that the bulk viscosity vanishes: 
( = 0. Hence, we shall concentrate our attention on the shear viscosity r\. 

2.3 Viscosity In Weakly Coupled Field Theories 

We now briefly consider the behavior of the shear viscosity in weakly coupled field theories, with 
the X(j) 4 theory as a concrete example. At weak coupling, there is a separation between two length 
scales: The mean free path of particles is much larger than the distance scales over which scatterings 
occur. Each scattering event takes a time of order T -1 (which can be thought of as the time required 
for final particles to become on-shell) . The mean free path ^ m f p can be estimated from the formula 

Cf P ~ — , (24) 
nav 

where n is the density of particles, a is the typical scattering cross section, and v is the typical 
particle velocity. Inserting the values for thermal A0 4 theory, n ~ T 3 , a ~ \ 2 T~ 2 , and v ~ 1, one 
finds . 

4nfp ~ ^ » T ' (25) 

The viscosity can be estimated from kinetic theory to be 



(26) 
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where e is the energy density. Prom e ~ T and the estimate of £ m ft, one finds 

T 3 

V~#- (27) 

In particular, the weaker the coupling A, the larger the viscosity n. This behavior is explained by 
the fact that the viscosity measures the rate of momentum diffusion. The smaller A is, the longer 
a particle travels before colliding with another one, and the easier the momentum transfer. 

It may appear counterintuitive that viscosity tends to infinity in the limit of zero coupling A — > 0: 
At zero coupling there is no dissipation, so should the viscosity be zero? The confusion arises owing 
to the fact that the hydrodynamic theory, and hence the notion of viscosity, makes sense only on 
distances much larger than the mean free path of particles. If one takes A —* 0, then to measure 
the viscosity one has to do the experiment at larger and larger length scales. If one fixes the size 
of the experiment and takes A — > 0, dissipation disappears, but it does not tell us anything about 
the viscosity. 

As will become apparent below, a particularly interesting ratio to consider is the ratio of shear 
viscosity and entropy density s. The latter is proportional to T 3 ; thus 

<*) 

s \ z 

One has i//s>1 for A<1. This is a common feature of weakly coupled field theories. Extrapolat- 
ing to A ~ 1, one finds r\j s ~ 1. We shall see that theories with gravity duals are strongly coupled, 
and rj/s is of order one. More surprisingly, this ratio is the same for all theories with gravity duals. 

To compute rather than estimate the viscosity, one can use Kubo's formula. It turns out that 
one has to sum an infinite number of Feynman graphs to even find the viscosity to leading order. 
Another way that leads to the same result is to first formulate a kinetic Boltzmann equation for the 
quasi-particles as an intermediate effective description, and then derive hydrodynamics by taking 
the limit of very long lengths and time scales in the kinetic equation. Interested readers should 
consult Refs. 115, 161 for more details. 



3 AdS/CFT CORRESPONDENCE 

3.1 Review Of AdS/CFT Correspondence At Zero Temperature 

This section briefly reviews the AdS/CFT correspondence at zero temperature. It contains only the 
minimal amount of materials required to understand the rest of the review. Further information 
can be found in existing reviews and lecture notes [T7][T5]. 

The original example of AdS/CFT correspondence is between N = 4 supersymmetric Yang-Mills 
(SYM) theory and type IIB string theory on AdSsxS 5 space. Let us describe the two sides of the 
correspondence in some more detail. 

The M = 4 SYM theory is a gauge theory with a gauge field, four Weyl fermions, and six real 
scalars, all in the adjoint representation of the color group. Its Lagrangian can be written down 
explicitly, but is not very important for our purposes. It has a vanishing beta function and is a 
conformal field theory (CFT) (thus the CFT in AdS/CFT). In our further discussion, we frequently 
use the generic terms "field theory" or CFT for the M = 4 SYM theory. 
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On the string theory side, we have type IIB string theory, which contains a finite number of 
massless fields, including the graviton, the dilaton <5, some other fields (forms) and their fermionic 
superpartners, and an infinite number of massive string excitations. It has two parameters: the 
string length l s (related to the slope parameter a' by a' = I 2 ) and the string coupling g s . In the 
long- wavelength limit, when all fields vary over length scales much larger than l s , the massive modes 
decouple and one is left with type IIB supergravity in 10 dimensions, which can be described by 
an action [19J 

Ssugra = ttV / dl ° x ^9 e~ 2 * {K + 4 d^d^ + •••), (29) 

where K\o is the 10-dimensional gravitational constant, 

Kio = v / 8^G = 8vr 7/2 5s ^ , (30) 

and • • • stay for the contributions from fields other than the metric and the dilaton. One of these 
fields is the five- form F§, which is constrained to be self-dual. The type IIB string theory lives is a 
10-dimensional space-time with the following metric: 

r 2 R 2 
ds 2 = -^(-dt 2 + dx. 2 ) + ^dr 2 + R 2 dn 2 5 . (31) 

The metric is a direct product of a five-dimensional sphere (dQ, 2 ) and another five-dimensional 
space-time spanned by t, x, and r. An alternative form of the metric is obtained from Eq. (|3ip by 
a change of variable z = R 2 /r, 



r>2 

ds 2 = — {-dt 2 + dx 2 + dz 2 )+R 2 dnl. (32) 



Both coordinates r and z are known as the radial coordinate. The limiting value r = oo (or z = 0) 
is the boundary of the AdS space. 

It is a simple exercise to check that the (i, x, r) part of the metric is a space with constant 
negative curvature, or an anti de-Sitter (AdS) space. To support the metric (f3Tj) (i.e., to satisfy the 
Einstein equation) there must be some background matter field that gives a stress-energy tensor in 
the form of a negative cosmological constant in AdSs and a positive one in S 5 . Such a field is the 
self-dual five-form field F§ mentioned above. 

Field theory has two parameters: the number of colors N and the gauge coupling g. When the 
number of colors is large, it is the 't Hooft coupling A = g 2 N that controls the perturbation theory. 
On the string theory side, the parameters are g s , l s , and radius R of the AdS space. String theory 
and field theory each have two dimensionless parameters which map to each other through the 
following relations: 

9 2 = 4tt 5s , (33) 
g 2 N c = (34) 

Equation (|33p tells us that, if one wants to keep string theory weakly interacting, then the gauge 
coupling in field theory must be small. Equation (|34p is particularly interesting. It says that the 
large 't Hooft coupling limit in field theory corresponds to the limit when the curvature radius of 
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space-time is much larger than the string length l s . In this limit, one can reliably decouple the 
massive string modes and reduce string theory to supergravity. In the limit g s <C 1 and R S> l s , 
one has classical supergravity instead of string theory. The practical utility of the AdS/CFT 
correspondence comes, in large part, from its ability to deal with the strong coupling limit in gauge 
theory. 

One can perform a Kaluza-Klein reduction [20] by expanding all fields in S 5 harmonics. Keeping 
only the lowest harmonics, one finds a five-dimensional theory with the massless dilaton, SO(6) 
gauge bosons, and gravitons [2T] : 

N 2 f ( 1 R 2 \ 

Ssb = ^3 J d " x [ n ™ - 2A - -d^d^ - —F^F^ + •••). (35) 

In AdS/CFT, an operator O of field theory is put in a correspondence with a field cj) ("bulk" field) 
in supergravity. We elaborate on this correspondence below; here we keep the operator and the field 
unspecified. In the supergravity approximation, the mathematical statement of the correspondence 
is 

Z m [J] = e iS[M . (36) 

On the left is the partition function of a field theory, where the source J coupled to the operator 
O is included: 

Z AD [J] = lD4> exp (iS + iJ d 4 x JO) . (37) 

On the right, S'f^d] is the classical action of the classical solution <j) c \ to the field equation with the 
boundary condition: 

lim ^ = J{x) . (38 ) 

z^O Z A 

Here A is a constant that depends on the nature of the operator O (namely, on its spin and 
dimension). In the simplest case, A = 0, and the boundary condition becomes (j) c \{z=Q) = J. 
Differentiating Eq. (|36p with respect to J, one can find the correlation functions of O. For example, 
the two-point Green's function of O is obtained by differentiating S c i[</>] twice with respect to the 
boundary value of 4>, 



G(x - y) = -i(TO(x)0(y)) = ^ S ' [<V|] 



SJ(x)SJ(y) 



(39) 

=0)=J 



The AdS/CFT correspondence thus maps the problem of finding quantum correlation functions in 
field theory to a classical problem in gravity. Moreover, to find two-point correlation functions in 
field theory, one can be limited to the quadratic part of the classical action on the gravity side. 

The complete operator to field mapping can be found in Refs. [5JH7]. For our purpose, the 
following is sufficient: 

• The dilaton $ corresponds to O = —C = \F 2 U H , where C is the Lagrangian density. 

• The gauge field corresponds to the conserved R-charge current J a ^ of field theory. 

• The metric tensor corresponds to the stress-energy tensor T^ v . More precisely, the partition 
function of the four-dimensional field theory in an external metric g® is equal to 

^4d[s°J = exp(iS d [?/«/]) , (40) 
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where the five-dimensional metric satisfies the Einstein's equations and has the following 
asymptotics at z = 0: 

ds 2 = g^dx^dx" = -^{dz 2 + g^dx^dx") . (41) 

From the point of view of hydrodynamics, the operator jF 2 is not very interesting because its 
correlator does not have a hydrodynamic pole. In contrast, we find the correlators of the R-charge 
current and the stress-energy tensor to contain hydrodynamic information. 

We simplify the graviton part of the action further. Our two-point functions are functions of 
the momentum p = (u>,k). We can choose spatial coordinates so that k points along the x 3 -axis. 
This corresponds to perturbations that propagate along the x 3 direction: h^ u = h^ii, r, x 3 ). These 
perturbations can be classified according to the representations of the 0(2) symmetry of the (x , x 2 ) 
plane. Owing to that symmetry, only certain components can mix; for example, h\2 does not mix 
with any other components, whereas components hoi and /131 mix only with each other. We assume 
that only these three metric components are nonzero and introduce shorthand notations 

<t> = h\, oo = hi, a 3 = h\. (42) 

The quadratic part of the graviton action acquires a very simple form in terms of these fields: 

N 2 r ( 1 1 \ 

S qU ad = ^3 Jd 4 xdr^(--g^d^d^-^-g^g^f^f a A , (43) 

where = d^a u — d u a^, and g 2 s = g xx . In deriving Eq. (filjj) . our only assumption about the 
metric is that it has a diagonal form, 

ds 2 = g tt dt 2 + g„dr 2 + g xx dx 2 , (44) 

so it can also be used below for the finite-temperature metric. 

As a simple example, let us compute the two-point correlation function of T xy , which corresponds 
to (j) in gravity. The field equation for (j) is 

d^V=g~g^d u cf ) ) = 0. (45) 
The solution to this equation, with the boundary condition <j)(p, z = 0) = 4>o(p), can be written as 

<f>(p,z) = f p (z)<h(p), (46) 
where the mode function f p (z) satisfies the equation 



with the boundary condition / p (0) = 1. The mode equation (|47p can be solved exactly. Assuming 
p is spacelike, p 2 > 0, the exact solution and its expansion around z = is 

f p (z) = \{pz) 2 K 2 {pz) = 1 - -ipz) 2 - UpzfHpz) + 0{{pzf) . (48) 
2 4 lb 

The second solution to Eq. (|4T|) . (pz) 2 l2(pz), is ruled out because it blows up at z — > 00. 
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We now substitute the solution into the quadratic action. Using the field equation, one can 
perform integration by parts and write the action as a boundary integral at z = 0. One finds 

N 2 f , 1 , f d 4 p 

S = Wtt 2 J d x ^3^( x ' z )^( x ' z )l^o = J 7^4 ^o(-p)^ r (p,2r)0o(p)U-*o J (49) 

where 2 

Hp, z) = ^2^f-P^)9zf P (z) ■ (50) 
Differentiating the action twice with respect to the boundary value 4>q one finds 

AT 2 

(T xy T xy ) p = -21imF(p,z) = — p 4 ln(p 2 ) . (51) 

Note that we have dropped the term ~ p 4 lnz, which, although singular in the limit z — ► 0, is a 
contact term [i.e., a term proportional to a derivative of 5{x) after Fourier transform]. Removing 
such terms by adding local counter terms to the supergravity action is known as the holographic 
renormalization [22J. It is, in a sense, a holographic counterpart to the standard renormalization 
procedure in quantum field theory, here applied to composite operators. 

For time-like p, p 2 < 0, there are two solutions to Eq. (|47p which involve Hankel functions H^\z) 
and H( 2 \z) instead of K2(z). Neither function blows up at z — > oo, and it is not clear which should 
be picked. Here we encounter, for the first time, a subtlety of Minkowski-space AdS/CFT, which is 
discussed in great length in subsequent sections. At zero temperature this problem can be overcome 
by an analytic continuation from space-like p. However, this will not work at nonzero temperatures. 

3.2 Black Three-Brane Metric 

At nonzero temperatures, the metric dual to M = 4 SYM theory is the black three-brane metric, 

ds 2 = ^(-fdt 2 + (fx 2 ) + ^Ur 2 + R 2 dn 2 , (52) 
Rr r z j 

with / = 1 — rg/r 4 . The event horizon is located at r = ro, where / = 0. In contrast to the usual 
Schwarzschild black hole, the horizon has three flat directions x. The metric (I52p is thus called a 
black three-brane metric. 

We frequently use an alternative radial coordinate u, defined as u = r^/r 2 . In terms of u, the 
boundary is at u = 0, the horizon at u = 1, and the metric is 

dg2 = {KTRf 2 + ^ + _& du 2 + R 2 m 2 (53) 

u z Au z j{u) 

The Hawking temperature is determined completely by the behavior of the metric near the 
horizon. Let us concentrate on the (t, r) part of the metric, 

R A 4ro(r — ro) 

Changing the radial variable from r to p, 

P 2 

r = r H , (55) 
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and the metric components become nonsingular: 

ds 2 = ^ [dp 2 - ^P 2 dt 2 ^ . (56) 

Note also that after a Wick rotation to Euclidean time r, the metric has the form of the flat metric 
in cylindrical coordinates, ds 2 ~ dp 2 + p 2 d<p 2 , where (p = 2tqR~ 2 t. To avoid a conical singularity 
at p = 0, ip must be a periodic variable with periodicity 2ir. This fact matches with the periodicity 
of the Euclidean time in thermal field theory r ~ r + 1/T, from which one finds the Hawking 
temperature: 

Th = ^- (57) 

One of the first finite-temperature predictions of AdS/CFT correspondence is that of the ther- 
modynamic potentials of the N = 4 SYM theory in the strong coupling regime. The entropy is 
given by the Bekenstein-Hawking formula S = A /{AG), where A is the area of the horizon of the 
metric ([52]) ; the result can then be converted to parameters of the gauge theory using Eqs. ([30]) . 
([M]), and (J31D. One obtains 

S n2 -N 2 T 3 , (58) 



V 2 

which is 3/4 of the entropy density in M = 4 SYM theory at zero 't Hooft coupling. 

We now try to generalize the AdS/CFT prescription to finite temperature. In the Euclidean 
formulation of finite-temperature field theory, field theory lives in a space-time with the Euclidean 
time direction r compactified. The metric is regular at r = tq: If one views the (r, r) space as 
a cigar-shaped surface, then the horizon r = r$ is the tip of the cigar. Thus, ro is the minimal 
radius where the space ends, and there is no point in space with r less than ro. The only boundary 
condition at r = ro is that fields are regular at the tip of the cigar, and the AdS/CFT correspondence 
is formulated as 

Z 4 d[J] = ^5d[0]|^= O )^j • (59) 



4 REAL-TIME AdS/CFT 

In many cases we must find real-time correlation functions not given directly by the Euclidean path- 
integral formulation of thermal field theory. One example is the set of kinetic coefficients expressed, 
through Kubo's formulas, via a certain limit of real-time thermal Green's functions. Another related 
example appears if we want to directly find the position of the poles in the correlation functions 
that would correspond to the hydrodynamic modes. 

In principle, some real-time Green's functions can be obtained by analytic continuation of the 
Euclidean ones. For example, an analytic continuation of a two-point Euclidean propagator gives 
a retarded or advanced Green's function, depending on the way one performs the continuation. 
However, it is often very difficult to directly compute a quantity of interest in that way. In particular, 
it is very difficult to get the information about the hydrodynamic (small u>, small k) limit of real-time 
correlators from Euclidean propagators. The problem here is that we need to perform an analytic 
continuation from a discrete set of points in Euclidean frequencies (the Matsubara frequencies) 
u = 2irin, where n is an integer, to the real values of to. In the hydrodynamic limit, we are 
interested in real and small u, whereas the smallest Matsubara frequency is already 2ttT. 
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Therefore, we need a real-time AdS/CFT prescription that would allow us to directly compute 
the real-time correlators. However, if one tries to naively generalize the AdS/CFT prescription, one 
immediately faces a problem. Namely, now r = r$ is not the end of space but just the location of 
the horizon. Without specifying a boundary condition at r = ro, there is an ambiguity in defining 
the solution to the field equations, even as the boundary condition at r = oo is set. 

As an example, let us consider the equation of motion of a scalar field in the black hole back- 
ground, d^{g^ v d u <p) = 0. The solution to this equation with the boundary condition <f> = 0o at 
u = is (j)(p, u) = f p (f>o(p), where f p (u) satisfies the following equation in the metric ([53]) : 



1 + u 2 „, w 2 



•2 



fl-^fp + ^fp- -JP = *- (60) 

Here the prime denotes differentiation with respect to u, and we have defined the dimensionless 
frequency and momentum: 

w = — ~ i Q= — — ■ 61) 
2?rT h 2ttT V ; 

Near u = the equation has two solutions, /i ~ 1 and f'2 ~ u 2 . In the Euclidean version of 

thermal AdS/CFT, there is only one regular solution at the horizon u = 1, which corresponds to a 

particular linear combination of f± and f%. However, in Minkowski space there are two solutions, 

and both are finite near the horizon. One solution termed f p behaves as (1 — u)~ lw / 2 , and the other 

is its complex conjugate /* ~ (1 — u) tw ^ 2 . These two solutions oscillate rapidly as u — > 1, but the 

amplitude of the oscillations is constant. Thus, the requirement of finiteness of f p allows for any 

linear combination of /1 and f% near the boundary, which means that there is no unique solution 

to Eq. (1601) . 



where the coordinate 



4.1 Prescription For Retarded Two-Point Functions 

Physically, the two solutions f p and /* have very different behavior. Restoring the e~' lujt phase in 
the wave function, one can write 

e'^fp ~ e-' 1 ^^ , (62) 
e~ lujt f; ~ e ^(*- r *) , (63) 

ln(l — u) , 

was introduced so that Eqs. (|62|) and (|63|) looked like plane waves. In fact, Eq. (|62|) corresponds to 
a wave that moves toward the horizon (incoming wave) and Eq. (|63p to a wave that moves away 
from the horizon (outgoing wave). 

The simplest idea, which is motivated by the fact that nothing should come out of a horizon, 
is to impose the incoming-wave boundary condition at r = rg and then proceed as instructed by 
the AdS/CFT correspondence. However, now we encounter another problem. If we write down 
the classical action for the bulk field, after integrating by parts we get contributions from both the 
boundary and the horizon: 

d^p z=z h 
■ <t)n(— v^Fiv, z)d>n(v) _ q 



s = I ^-j4 M-p)f(p> z )<h(p) . „ • ( 65 ) 
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If one tried to differentiate the action with respect to the boundary value 4>q, one would find 

G(p)=r(p,z)\ z B +F(-p,z)\ s H . (66) 

From the equation satisfied by f p and from /* = f_ p , it is easy to show that the imaginary part 
of J-(p, z) does not depend on z\ hence the quantity G(p) in Eq. (|66|) is real. This is clearly not 
what we want, as the retarded Green's functions are, in general, complex. Simply throwing away 
the contribution from the horizon does not help because 0~{— p, z) = J-*(p,z) owing to the reality 
of the equation satisfied by f p . 

A partial solution to this problem was suggested in Ref. p 7 ]. It was postulated that the re- 
tarded Green's function is related to the function T by the same formula that was found at zero 
temperature: 

G R (p) = -21imf(p,z). (67) 

In particular, we throw away all contributions from the horizon. This prescription was established 
more rigorously in Ref. [5] (following an earlier suggestion in Ref. [23]) as a particular case of a 
general real-time AdS/CFT formulation, which establishes the connection between the close-time- 
path formulation of real-time quantum field theory with the dynamics of fields in the whole Penrose 
diagram of the AdS black brane. Here we accept Eq. (|57|) as a postulate and proceed to extract 
physical results from it. 

It is also easy to generalize this prescription to the case when we have more than one field. In 
that case, the quantity T becomes a matrix J 7 ^, whose elements are proportional to the retarded 
Green's function G^. 



4.2 Calculating Hydrodynamic Quantities 

As an illustration of the real-time AdS/CFT correspondence, we compute the correlator of T xy . 
First we write down the equation of motion for <f> = hyi 

„ 1 + u 2 , w 2 - q 2 f 
<\> p ~ -^j-^v + u p = • ( 68 ) 

In contrast to the zero-temperature equation, now u and k enter the equation separately rather 
than through the combination oj 2 — k 2 . Thus the Green's function will have no Lorentz invariance. 
The equation cannot be solved exactly for all u and k. However, when u and k are both much 
smaller than T, or w,q <C 1, one can develop series expansion in powers of w and q. There are 
two solutions that are complex conjugates of each other. The solution that is an incoming wave at 
u = 1 and normalized to 1 at u = is 

f p (z) = (l-u 2 )-™/ 2 + 0(w 2 ,q 2 ). (69) 

The kinetic term in the action for 4> is 

S—^(*fj>. (70) 



u 



Applying the general formula (|67p . one finds the retarded Green's function of T ; 



G*{u, k) = - " - A - iw , (71) 



ir 2 N 2 T A 



■xy,xy 
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and, using Kubo's formula for rj, the viscosity, 



7T 



r, = —N T 6 . (72) 
8 

It is instructive to compute other correlators that have poles corresponding to hydrodynamic 
modes. As a warm-up, let us compute the two-point correlators of the R-charge currents, which 
should have a pole at u> = —iDk 2 , where D is the diffusion constant. We first write down Maxwell's 
equations for the bulk gauge field. Let the spatial momentum be aligned along the x 3 -axis: p = 
(to, 0, 0, k). Then the equations for Aq and ^3 are coupled: 

wA' + qfA' 3 = 0, (73) 
A'^-—(q 2 A Q + wqA 3 ) = 0, (74) 



uf 

.\" + jA^ + ^(w 2 A 3 ^,r ( ,.\ () ) = 0. (75) 



One can eliminate ^3 and write down a third-order equation for Aq, 

K + ^A>> + ^-^A> = 0. (76) 

Near u=l we find two independent solutions, A' Q ~ (1 — u)^ lw l 2 , and the incoming-wave boundary 
condition singles out (1 - u)' iw/2 . One can substitute A' = (1 - u)~ iw / 2 F(u) into Eq. (JT6D . The 
resulting equation can be solved perturbatively in w and q 2 . We find 



(78) 



Using Eq. (1741) one can express C through the boundary values of Aq and ^3 at u = 0: 

q 2 A + wqA 3 



C 



iw — q 2 



u=0 

Differentiating the action with respect to the boundary values, we find, in particular, 

N 2 T k 2 

Wo) P = -r r -. p^, (79) 

167r iuj — Dk z 

where 



The correlator given by Eq. ([79j) has the expected hydrodynamic diffusive pole, and D is the 
R-charge diffusion constant. 

Similarly, one can observe the appearance of the shear mode in the correlators of the metric 
tensor. We note that the shear flow along the x 1 direction with velocity gradient along the x 3 
direction involves Tqi and T31, hence the interesting metric components are ao = h}) and 03 = h\. 
Two of the field equations are 

a' - q J-a' 3 = 0, (81) 
1 + u 2 1 

a'i —03 + —p; (w 2 a 3 + wqa ) = . (82) 

uf uf z 
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They can be combined into a single equation: 

2uf-q 2 f + w 2 , 



in 2m v 



«o 



a + 



/ ~ U ' up 
Again, the solution can be found perturbatively in w and q: 



a = 0. 



a = C(l - u)- iw/2 



u — 1 — it in — - — 

2 2 



Applying the prescription, one finds the retarded Green's functions. For example, 

ik 2 



where 



7T 



£ = -N 2 T 3 , 



iuj — Vk 2 
1 



V 



(83) 



(84) 



(85) 



(86) 



Thus, we found that the correlator contains a diffusive pole u = —iDk 2 , just as anticipated from 
hydrodynamics. Furthermore, the magnitude of the momentum diffusion constant T> also matched 
our expectation. Indeed, if one recalls the value of rj from Eq. (|72p and the entropy density from 
Eq. ([58]) . one can check that 

V = — ^ . (87) 



e + P 



5 THE MEMBRANE PARADIGM 



Let us now look at the problem from a different perspective. The existence of hydrodynamic modes 
in thermal field theory is reflected by the existence of the poles of the retarded correlators computed 
from gravity. Are there direct gravity counterparts of the hydrodynamic normal modes? 

If the answer to this question is yes, then there must exist linear gravitational perturbations of 
the metric that have the dispersion relation identical to that of the shear hydrodynamic mode, 
uj ~ — iq 2 , and of the sound mode, uj = c s q — i^q 2 . It turns out that one can explicitly construct the 
gravitational counterpart of the shear mode. (It should be possible to find a similar construction 
for the sound mode, but it has not been done in the literature; for a recent work on the subject, 
see [21] •) Our discussion is physical but somewhat sketchy; for more details see Ref. [25] . 

First, let us construct a gravity perturbation that corresponds to a diffusion of a conserved charge 
(e.g., the R-charge in M = 4 SYM theory). To keep the discussion general, we use the form of the 
metric ()44p . with the metric components unspecified. Our only assumptions are that the metric is 
diagonal and has a horizon at r = ro, near which 

goo = -7o 0* - ?~o), 9rr = — — • (88) 

r — ro 

The Hawking temperature can be computed by the method used to arrive at Eq. (|5T|) . and one 
finds T = (4^)- 1 (7o/7r) 1/2 - 

We also assume that the action of the gauge field dual to the conserved current is 

Sgauge = Jdx^g- ^-^-F^F^ , (89) 



Viscosity, Black Holes, and QFT 



17 



where g e R is an effective gauge coupling that can be a function of the radial coordinate r. For 
simplicity we set g e R to a constant in our derivation of the formula for D; it can be restored by 
replacing y/—g — > \/ — ff/fleff m the final answer. 
The field equations are 

1 




-gF^\ = . (90) 

We search for a solution to this equation that vanishes at the boundary and satisfies the incoming- 
wave boundary condition at the horizon. 

The first indication that one can have a hydrodynamic behavior on the gravity side is that 
Eq. (|90p implies a conservation law on a four-dimensional surface. We define the stretched horizon 
as a surface with constant r just outside the horizon, 

r = r h = r + e, e <C r , (91) 

and the normal vector directed along the r direction (i.e., perpendicularly to the stretched 
horizon). Then with any solution to Eq. (|90p . one can associate a current on the stretched horizon: 

f = n v F^ . (92) 

The antisymmetry of F^ v implies that j^ has no radial component, j r = 0. The field equation (|90j) 
and the constancy of n v on the stretched horizon imply that this current is conserved: duj^ = 0. To 
establish the diffusive nature of the solution, we must show the validity of the constitutive equation 

ji = -Ddd°. 

Such constitutive equation breaks time reversal and obviously must come from the absorptive 
boundary condition on the horizon. The situation is analogous to the propagation of plane waves to 
a non-reflecting surface in classical electrodynamics. In this case, we have the relation B = — n x E 
between electric and magnetic fields. In our case, the corresponding relation is 

Fir = (93) 
V 7o r - r 

valid when r is close to rp. This relates ji ~ F{ r to the parallel to the horizon component of the 
electric field Fqi, which is one of the main points of the "membrane paradigm" approach to black 
hole physics [26[127j- We have yet to relate ji to jo ~ F$ r , which is the component of the electric 
field normal to the horizon. To make the connection to F$ r , we use the radial gauge A r = 0, in 
which 

F 0l « -dtAo . (94) 

Moreover, when k is small the fields change very slowly along the horizon. Therefore, at each 
point on the horizon the radial dependence of the scalar potential Aq is determined by the Poisson 
equation, 

d r (V=gg rr g 00 d r A ) = o, (95) 

whose solution, which satisfies Aq(t = oo) = 0, is 

oo 

Mr) = cjdr> 9 ^; )9 ^ ) . (96) 

J v-g( r ) 
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This means that the ratio of the scalar potential Aq and electric field Fq t approaches a constant 
near the horizon: 

oo 

A 



F 



Or 



r=rg 



- g -{r Q ) /dr^fer). (97) 



9009rr J \/—g 



Combining the formulas f ~ F^ ~ O-iAq, and A$ ~ Fq t ~ j°, we find Fick's law f = —Ddij , with 
the diffusion constant 

oo 2 

r> V "5 / \ f r -9009rrg c s , \ /nn\ 

D = g — / ( r o) / rfr / — ( r ) • ( 98 ) 

r 

Thus, we found that for a slowly varying solution to Maxwell's equations, the corresponding 
charge on the stretched horizon evolves according to the diffusion equation. Therefore, the gravity 
solution must be an overdamped one, with ui = —iDk 2 . This is an example of a quasi-normal mode. 
We also found the diffusion constant D directly in terms of the metric and the gauge coupling g e g . 

The reader may notice that our quasinormal modes satisfy a vanishing Dirichlet condition at 
the boundary r=oo. This is different from the boundary condition one uses to find the retarded 
propagators in AdS/CFT, so the relation of the quasinormal modes to AdS/CFT correspondence 
may be not clear. It can be shown, however, that the quasi-normal frequencies coincide with the 
poles of the retarded correlators [28j l29]. 

We can now apply our general formulas to the case of M = 4 SYM theory. The metric components 
are given by Eq. (|52p . For the R-charge current g c R = const, Eq. (|98p gives D = l/(2irT), in 
agreement with our AdS/CFT computation. For the shear mode of the stress-energy tensor we 
have effectively g\ s = g XX: so V = l/(47rT), which also coincides with our previous result. In both 
cases, the computation is much simpler than the AdS/CFT calculation. 



6 THE VISCOSITY/ENTROPY RATIO 
6.1 Universality 

In all thermal field theories in the regime described by gravity duals the ratio of shear viscosity r] 
to (volume) density of entropy s is a universal constant equal to l/(47r) [h/ '(47tA:b), if one restores 
h, c and the Boltzmann constant fee]. 

One proof of the universality is based on the relationship between graviton's absorption cross 
section and the imaginary part of the retarded Green's function for T xy [31]. Another way to prove 
the universality [32] is via the direct AdS/CFT calculation of the correlation function in Kubo's 
formula Iflg]) . 

We, however, follow a different method. It is based on the formula for the viscosity derived from 
the membrane paradigm. A similar proof was given by Buchel & Liu [30J . 

The observation is that the shear gravitational perturbation with k = can be found exactly by 
performing a Lorentz boost of the black-brane metric f)52[) . Consider the coordinate transformations 
r, t, Xi — ► r', t' , x[ of the form 



r 
t 
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vt' 



y' + vt' , 



(99) 



where v < 1 is a constant parameter and the expansion on the right corresponds to v <C 1. In the 
new coordinates, the metric becomes 



ds 2 = goodt' 2 + g„dr 2 + g xx {r) ^{dx l ) 2 + 2v(g 00 + g xx )dt' dy' 



(100) 



8=1 



This is simply a shear fluctuation at k = 0. In our language, the corresponding gauge potential is 



a 



vg xx (g Q0 + g x 



(101) 



This field satisfies the vanishing boundary condition ao(r = oo) = owing to the restoration 
of Poincare invariance at the boundary: goo/g xx — * — 1 when r — > oo. This clearly has a much 
simpler form than Eq. (196j) for the solution to the generic Poisson equation. The simple form of 



solution (jlOip is valid only for the specific case of the shear gravitational mode with g^ = g xx . 
We have also implicitly used the fact that the metric satisfies the Einstein equations, with the 
stress-energy tensor on the right being invariant under a Lorentz boost. 
Equation (|97|) now becomes 



The shear mode diffusion constant is 



V = -± 

JOr 



i + g xx goo 



d r (g xx goo) 



9xx(ro) 
7o 



■»■'■() 



V707r 

fc(ro) 



4vrT 



Because V = rj/(e + P), and e + P 



Ts in the absence of chemical potentials, we find that 
r/ 1 

S 47T 



(102) 



(103) 



(104) 



In fact, the constancy of this ratio has been checked directly for theories dual to Dp-brane 
M-brane [11], Klebanov-Tseytlin and Maldacena-Nunez backgrounds [30], N = 2* SYM theory 
[33] and others. Curiously, the viscosity to entropy ratio is also equal to 1/47T in the pre-AdS/CFT 
"membrane paradigm" hydrodynamics [M]: there, for a four-dimensional Schwarzschild black hole 
one has rj m _ p . = I/IGttGn, while the Bekenstein-Hawking entropy is s = 1/4Gjv. 

As remarked in Sec. [2l the ratio rj/s is much larger than the one for weakly coupled theories. 
The fact that we found the ratio to be parametrically of order one implies that all theories with 
gravity duals are strongly coupled. 

In J\f = 4 SYM theory, the ratio n/s has been computed to the next order in the inverse 't Hooft 
coupling expansion [3S] 

r, 1 / 135C(3) \ 

s = I 1 + wm 2 ) ■ ( } 

The sign of the correction can be guessed from the fact that in the limit of zero 't Hooft coupling 
g 2 N — > 0, the ratio diverges, rj/s — > oo. 
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From our discussion above, one can argue that 

2 > f (106) 

in all systems that can be obtained from a sensible relativistic quantum field theory by turning on 
temperatures and chemical potentials. 

The bound, if correct, implies that a liquid with a given volume density of entropy cannot be 
arbitrarily close to being a perfect fluid (which has zero viscosity). As such, it implies a lower bound 
on the viscosity of the QGP one may be creating at RHIC. Interestingly, some model calculations 
suggest that the viscosity at RHIC may be not too far away from the lower bound [36l[37] . 

One place where one may think that the bound should break down is superfluids. The ability 
of a superfluid to flow without dissipation in a channel is sometimes described as "zero viscosity" . 
However, within the Landau's two-fluid model, any superfluid has a measurable shear viscosity 
(together with three bulk viscosities). For superfluid helium, the shear viscosity has been measured 
in a torsion-pendulum experiment by Andronikashvili [38]. If one substitutes the experimental 
values, the ratio rj/s for helium remains larger than h/A-irks ~ 6.08 x 10 -13 Ks for all ranges of 
temperatures and pressures, by a factor of at least 8.8. 

As discussed in Sec. 12.31 the ratio rj/s is proportional to the ratio of the mean free path and the 
de Broglie wavelength of particles, 

V ^mfp 



(107) 

s A 

For the quasi-particle picture to be valid, the mean free path must be much larger than the de 
Broglie wavelength. Therefore, if the coupling is weak and the system can be described as a 
collection of quasi-particles, the ratio rj/s is larger than 1. 

We have found is that, within the N = 4 SYM theory and, more generally, theories with gravity 
duals, even in the limit of infinite coupling the ratio rj/s cannot be made smaller than l/(4-7r). 



7 CONCLUSION 

In this review, we covered only a small part of the applications of AdS/CFT correspondence to 
finite-temperature quantum field theory. Here we briefly mention further developments and refer 
the reader to the original literature for more details. 

In addition to N = 4 SYM theory, there exists a large number of other theories whose hydrody- 
namic behavior has been studied using the AdS/CFT correspondence, including the worldvolume 
theories on M2- and M5-branes [11], theories on Dp branes [25], and little string theory [39]. In 
all examples the ratio rj/s is equal to l/(47r), which is not surprising because the general proofs of 
Sec. [6] apply in these cases. 

We have concentrated on the shear hydrodynamic mode, which has a diffusive pole (uj ~ —ik 2 ). 
One can also compute correlators which have a sound-wave pole from the AdS/CFT prescription 
|10j . One such correlator is between the energy density T 00 at two different points in space-time. 
The result confirms the existence of such a pole, with both the real part and imaginary part having 
exactly the values predicted by hydrodynamics (recall that in conformal field theories the bulk 
viscosity is zero and the sound attenuation rate is determined completely by the shear viscosity). 
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Some of the theories listed above are conformal field theories, but many are not (e.g., the Dp- 
brane worldvolume theories with p ^ 3). The fact that 77/s = l/(4-7r) also in those theories implies 
that the constancy of this ratio is not a consequence of conformal symmetry. Theories with less 
than maximal number of supersymmetries have been found to have the universal value of rj/s, for 
example, the N = 2* theory [40J , theories described by Klebanov-Tseytlin, and Maldacena-Nunez 
backgrounds [30J. A common feature of these theories is that they all have a gravitational dual 
description. The bulk viscosity has been computed for some of these theories |4iy39j. 

Besides viscosity, one can also compute diffusion constants of conserved charges by using the 
AdS/CFT correspondence. Above we presented the computation of the R-charge diffusion constant 
in J\f = 4 SYM theory; for similar calculations in some other theories see Ref. [Iiy25j. 

Recently, the AdS/CFT correspondence was used to compute the energy loss rate of a quark in 
the fundamental representation moving in a finite-temperature plasma [42, 43 , 44, 45 j. This quantity 
is of importance to the phenomenon of "jet quenching" in heavy- ion collisions. 

So far, the only quantity that shows a universal behavior at the quantitative level, across all 
theories with gravitational duals, is the ratio of the shear viscosity and entropy density. Recently, 
it was found that this ratio remains constant even at nonzero chemical potentials [46,47,48,49,50]. 

What have we learned from the application of AdS/CFT correspondence to thermal field theory? 
Although, at least at this moment, we cannot use the AdS/CFT approach to study QCD directly, 
we have found quite interesting facts about strongly coupled field theories. We have also learned 
new facts about quasi-normal modes of black branes. However, we have also found a set of puzzles: 
Why is the ratio of the viscosity and entropy density constant in a wide class of theories? Is there 
a lower bound on this ratio for all quantum field theories? Can this be understood without any 
reference to gravity duals? With these open questions, we conclude this review. 
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